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Abstract 

Integrated damping mechanics for composite laminates with constrained interlaminar layers 
of polymer damping materials are developed. Discrete layer damping mechanics are 
presented for composite laminates with damping layers, in connection with a semi-analytical 
method for predicting the modal damping in simply-supported specialty composite plates. 
Several application cases are used to demonstrate the advantages of the method. Damping 
predictions for graphite/epoxy composite plates of various laminations demonstrate the 
potential for higher damping than geometrically equivalent aluminum plates. The effects of 
aspect ratio, damping layer thickness, and fiber volume ratio on static and dynamic 
characteristics of the composite plate are also investigated. 
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Introduction 


Damping is a significant dynamic parameter for vibration and sound control, dynamic 
stability, positioning accuracy, fatigue endurance, and impact resistance. Many current 
structural applications (for example large space structures, engine blades, high speed 
machinery) require light-weight and high dynamic performance. Therefore, candidate sources 
of passive damping should add minimal parasitic weight and be compatible with the 
structural configuration. 

Two potential damping sources satisfying the previous requirements are the constrained 
damping layer approach, and the damping capacity of composites. Constrained damping 
layers in isotropic metallic structures have been widely applied and investigated 1 . They 
provide high damping, but tend to increase the structural weight and offer limited means for 
damping tailoring. The inherent damping capacity of composite materials also seems 
promising. Although the damping of composite structures is not veiy high, it is significantly 
higher than most common metallic structures. Moreover, composites are the materials of 
preference in many cases, as they readily provide high specific stiffness and strength. More 
importantly, research on the damping mechanics of composite laminates 2-4 and structures 5,6 
has shown that composite damping is anisotropic, highly tailorable, and depends on an array 
of micromechanical, laminate, and structural parameters. It has been further demonstrated 
that optimal tailoring may significantly improve the damped dynamic performance of 
composite structures 7 . 

It seems likely that the combination of both approaches, that is, composite structures with 
interlaminar damping layers will offer the advantages of high damping, damping tailoring, 
good mechanical properties, and low weight addition. In addition, the interlaminar damping 
concept is highly compatible with the laminated configuration of composite structures and 
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their fabrication techniques. Contrary to isotropic materials, the variations in anisotropy and 
elastic properties of each composite ply typically induce high interlaminar stresses, hence, 
interlaminar layers may produce higher and also tailorable damping. However, other critical 
mechanical properties, such as stiffness and strength, are expected to be reduced. In view 
of the complexity of the problem, it is very desirable to develop integrated composite 
mechanics that will efficiently predict the damping of composite structures with interlaminar 
damping layers, together with other critical mechanical properties and dynamic 
characteristics. Limited analytical and experimental work has been reported in this direction, 
including the damping of specialty composite truss members in extensional deformation 8 ' 10 , 
and the damping analysis of simple composite beams using continuum finite element 
simulations 11,12 . 

This paper presents recent developments on integrated damping mechanics for composite 
laminates with interlaminar damping layers. The laminate mechanics are general, in that they 
incorporate the capability for arbitrary numbers of plies and damping layers, arbitrary fiber 
orientations, and arbitrary fiber volume ratios (FVRs). They also include the capacity to 
simulate the damping of extensional, bending, and coupled deformations. The theory 
assumes a general piecewise continuous displacement field through the thickness of the 
laminate. The damping mechanics for composite laminates with damping layers is also 
integrated with other composite mechanics theories 13,14 , in order to include the effects of the 
many composite parameters on the dynamic characteristics. 

Finally, a semi-analytical methodology for the prediction of the dynamic properties of simply- 
supported (SS) specialty thick composite plates with interlaminar damping layers is 
developed. Results regarding the modal damping values, natural frequencies and static 
deflections of Gr/Epoxy composite plates are presented. The effects of various cross-ply 
configurations, aspect ratio, damping layer thickness, and FVR on the damped dynamic 
characteristics are investigated. 
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Method 


This section first reviews the synthesis of the composite properties for each composite ply 
and interlaminar layer. Then the developed damping mechanics for composite laminates 
with interply layers (Fig la), and the dynamic analysis of specialty composite plates with 
interlaminar damping layers are presented. 

Composite Properties 

For a unidirectional composite loaded along the material axes, closed-form expressions have 
been developed for the synthesis of elastic and dissipative properties 13,14 . Five independent 
elastic parameters completely characterize the stiffness of a unidirectional composite 
(orthotropic but transversely isotropic material): in-plane longitudinal modulus E IU , 
transverse in-plane modulus E l22 , in-plane Poisson’s ratio v 112 , in-plane shear modulus G )12 , 
and interlaminar shear modulus G ]23 . Additionally, four independent damping coefficients 
characterize the damping capacity, that is, longitudinal specific damping capacity (SDC) 
(direction 11), transverse in-plane damping »|; 12 (direction 22), transverse through-the- 
thickness damping (direction 33), in-plane shear damping i)> 16 (direction 12), 

interlaminar shear damping iji l4 (direction 23), and interlaminar shear damping «|» J6 
(direction 13). The specific damping capacities (SDCs) and elastic properties are explicitly 
related to fiber/matrix properties (ref. 14), and the fiber volume ratio (FVR). 

The interlaminar layers are assumed isotropic. Two independent elastic properties (Young’s 
modulus and Poisson’s ratio) and two independent damping properties (normal and shear 
SDC’s) completely characterize their viscoelastic behavior. 

Off-Axis Composites 
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For the case of off-axis composites, ie. composites loaded at an angle 0, the following 
transformations provide the stiffness matrix [EJ and the damping capacity matrix [vji c ] 
respectively 13,14 . 


[iy= [Rr l [E,\[Rr T 

(i) 

ng= M’twr 1, 

(2) 


[E,] is the stiffness matrix of the composite in the material coordinate system, and [4»,] 
represents the on-axis composite damping. The transformation matrix [R] and the damping 
matrices [»|/,] and [4> c ] are shown in the Appendix. 

Each interlaminar layer is treated as an isotropic material, hence, the above eqs. are reduced 
to trivial cases. In this manner, and contrary to many reported studies involving simplifying 
assumptions, all independent elastic and dissipative properties of the composite plies and the 
interlaminar laminae are considered. 

Laminates with Interlaminar Damping Layers 

To model the damping of laminates with interlaminar damping layers, a discrete layer 
laminate damping theory (DLDT) incorporating a piecewise continuous displacement field 
through-the-thickness is developed. Compared to other laminate theories based on assumed 
continuous displacement fields, the DLDT provides the capacity for more accurate strain 
and dissipative energy calculations in each ply, and particularly in each damping layer. The 
kinematic assumptions for the laminate theory are schematically shown in Fig lb. Discrete 
layer theories were originally proposed by Grigolyuk and were subsequently generalized by 
other researchers (see ref. 15) for the more accurate calculation of stresses in thick 
laminates. However, with additional development the DLDT combines the potential for 
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accurate damping predictions in composite laminates with damping layers, while maintaining 
generality and elegance. The assumed displacement field has the form, 


“(W >0= U °(W) + m(w>0 

v(xy^,t) = v°(xy,t) + v(xyj,t) 
w(xy?,t)= w°(xy,t) 

where superscript o represents the uniform through-the-thickness midplane deflection. 
Assuming that the displacements are separable functions of z, as suggested by Reddy and 
co-workers 16 , the previous eqs. take the form: 


N 

w(w,0= u°(xy,t)+ Y, ui { x ^,t)F\z) 

1 

N 

v(w ,/) = v °(xy,t) + Y vi(xy,t)F>(z) 


(4) 


where, u* and v* are displacements, along the x and y directions respectively, preferably at 
the interfaces between composite plies or sublaminates (group of plies) and the interlaminar 
damping layers. F*(z) are interpolation functions. In this manner, the assumed in-plane 
displacement field is general, in that it may represent extensional, flexural, shear, and 
coupled deformations, as well as, the interlaminar shear strains through-the-thickness of the 
laminate. 


The laminate strains are directly derived from eqs. (4) 


€ d = 4 + E 4^) i =1,2,6 

H 

€ a= 4 + E 4^( z ) i=4 > 5 

j’i 


(5) 
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where, the midplane strains are, 


if* 

ii 

* 

II 

s? 

0 0 

0 c 

II 

e c5=^ 


and the generalized strains are, 




II 

< 

K) 

II 

— * ^ 

ii 


II 

v> 

4= ui 




the comma in the subscripts indicates differentiation. 


( 6 ) 


( 7 ) 


The dissipated strain energy per unit area of the laminate Aw L is, 




1/2 r f/Kjwgc* 


J-W2 


( 8 ) 


Combination of eqs. (8, 5-7) ultimately provides the dissipated strain energy per unit area, 
Ah-, - V2 * 2 [S']< ♦ £ £ (9) 

y*i y-i i 

where, [A d ] is the extensional laminate damping matrix, including out-of-plane shear terms. 
The generalized coupling damping matrices [B d j ] and flexural/shear matrices [D d jm ] are new. 
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i/i =1,2,6 


N, 

Kl= E l- PM ’c¥z 

k- 1 J V/ 

N, 

(*k- E 

fe-i j **-j 
"/ * 

(B'X= E [ S >>=4,5 

k = l 

"t j ( 10 ) 

0>5V E [ - i/i =1,2,6 

k - 1 

N i 

<pf)*- e [- ^>= 4,5 

jt-i j **-; 


Similarly the maximum laminate strain energy per unit area w L is by definition, 


w 


L= 1/2 r «/[£„]«,* 


J- 7 i /2 


(ii) 


Combination of eqs. (11, 5-7) provides the maximum strain energy as a separable form of 
material properties and strains, v •• 

*> L - mCwi * ^fE mi * E E «?p*]0 ( 12 ) 

;'-l ;'-l ro-l 

where, [A] is the extensional laminate stiffness matrix, with additional out-of-plane shear 
terms. The generalized coupling stiffness matrices [B*] and flexural/shear matrices [iy m ] are 
new. 
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N, 

M- EL* 

*-l JA *-t 




E [** ([ E c \)^(z)dz 

k- 1 ih k-l 

iji =1,2,6 


N l 

c B %= E [ - m k )^)dz 

Ic. 1 Jh k-1 

i/t =4,5 

(13) 

N l 


(^)»- E f* 

k- 1 J *‘-' 

i/t =1,2,6 


N l 

( Djm )in= E [- 

*«1 Jh k-1 

i/i =4,5 



The equivalent laminate specific damping capacity for a given local displacement field is: 


i\> L = A wJw L 


(14) 


In both eqs. (10) and (13), the interlaminar damping layers are considered as individual plies 
having isotropic elastic and damping matrices [E c ] and [»Jj c ] respectively. For a given 
interlaminar layer, [4< c ] is determined from the mechanical properties of the material as a 
trivial composite case (ref. [14]). In this manner, the proposed laminate damping mechanics 
are general, as they can handle any number of composite plies, damping layers, and laminate 
configurations. 
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Simply-Supported Composite Plates 


The damped dynamic characteristics of simply supported specialty composite plates with 
interlaminar damping layers are calculated based on the previously described discrete-layer 
laminate damping theory. As will be shown in the remaining section, in the special case of 
composite plates consisting of 0 and 90 degrees plies, exact semi-analytical modal analysis 
solutions exist. Hence, exact and computationally inexpensive predictions of static and 
dynamic characteristics (modal damping, natural frequencies) are possible which will provide 
valuable insight in the mechanics of the problem and will establish the critical parameters. 
Most other laminations, boundary conditions, and structural configurations require 
approximate solutions, and such methods will be presented in the near future. 


Assuming a rectangular a by /? symmetric composite plate with negligible coupling (A 16 = 
A 26 =0, R> 16 = B i 26 =0, = D in, 26 =0), the following Navier fundamental solutions form 

a complete set of mode shapes in the x-y plane: 


*C.(W) = U° m cos(ax)sin(by)e ial 
= V L si n(ax)cos(by)e iut 

»*C.(W) = (b y y vl 

«i»(W) = U J mn cos(ax)sm(by)e ial 

v i»(w) = vL sin ( ax ) cos ( h y) eiwl 


(15) 


Where, a= mn/a and b= mr/fi. Combination of eqs. (6,7,15) yields the modal mid-plane and 
generalized strains as separable functions of x, y coordinates, time, and amplitudes. 

where, the terms in matrices [B] mn are sine/cosine functions of x, y coordinates and mode 
order. The amplitude displacement vectors are U° = {U°, V°, W°} T and U* = {U^, W} T . The 
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kinetic energy through the thickness of the plate is: 

K l = 1/2 j h ^{u} T p {a}dz (17) 

Considering eqs. (4) the laminate kinetic energy per unit area takes the form: 

k l - 1/2 * 2{a “) 7 i; 

J-1 ;'-l m-1 

where the generalized laminate mass and inertia matrices are 
"I . 

1 A m\ = E [ - dia 8(Pk) dz 

k-i 

N i 

[*il ' E L* diag(p„)Fl(z)dz 

*>i (19) 

N l 

[of 1 = E [*- 

The term diag(p k ) indicates a diagonal matrix, with all diagonal terms equal to the density 
of the k-th ply. 

By substituting eqs. (16) into the strain and kinetic energy expressions for the laminate, eqs. 
(12) and (18) respectively, integrating over the plate area, and applying Lagrangian 
dynamics, it can be proved that the undamped modal analysis (free vibration) solution of the 
plate takes the form: 

-to 2 \M 1 U + \K ]U =0 ( 20 ) 

where, the through-the-thickness modal displacements are {U mn } = {U omn ; U lmn ,...,U Nnin } 
and subscripts m, n indicate the mode order. Numerical solution of this eigenvalue problem 
provides the natural frequencies to mn and the through-the-thickness modes {U mn } for each 
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order ran of plane modes, in the context of eqs. (15). 


The modal SDC associated with the n-th vibration mode i|i n is: 


- 


J. 


A wdV 


J> 


( 21 ) 


dV 


where Aw n and w n represent the dissipated and maximum specific modal strain energies per 
cycle respectively. For laminated composite plates the modal damping becomes, where Aw Ln 


*. - 






( 22 ) 


and w Ln are respectively the dissipated and maximum laminate modal strain energies of the 
n-th mode. To calculate the modal damping, eqs. (10) and (13) are integrated through-the- 
thickness, and are combined with eqs. (9) and (12) respectively. Then the generalized 
stresses and strains are related to the calculated mode shape using eq. (16). The integrations 
over the area of the plate are performed analytically. Since the natural frequencies and the 
mode shapes through-the-thickness are calculated numerically from eq. (20), the whole 
procedure is semi-analytical. 
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Applications and Discussion 


This section presents applications of the developed methodology on square laminated 
Gr/epoxy simply-supported (SS) plates with interlaminar damping layers. The dimensions of 
the plate were assumed 11.8 in by 11.8 in. The composite material was HM-S graphite/epoxy 
with elastic and damping properties given in ref. 6. Unless otherwise stated, the fiber volume 
ratio (FVR) of the composite plies will be 0.50. The shear modulus and shear SDC of the 
interlaminar damping layer were 0.9 kpsi and 6.4 respectively, and represent measured 
values of a typical commercial damping polymer 17 at room temperature. Experimental data 
regarding the normal SDC of the damping polymer was not available, in part because their 
contributions are considered negligible. Hence, a low normal SDC (0.10) was assumed. The 
Poisson’s ratio of the damping polymer was assumed equal to 0.3. The effects of frequency 
on the properties of the damping polymer were neglected. The nominal thickness t f of the 
damping layer was equal to the typical thickness of a composite ply t, = 0.005 in. In all cases 
a concentrated vertical load of 100 lbs was applied at the center of the plate. Linear 
interpolation functions F j (z) were used to calculate the in-plane displacements through-the- 
thickness of the laminate. 

The ply-stacking sequence is represented using the traditional convention, but an 
interlaminar damping layer of thickness t ; is identified with symbol i. For comparison 
purposes, calculated results for aluminum plates are also presented, in such case, an 
aluminum layer of thickness equal to the thickness of k plies will be identified as Al k . All 
laminate configurations were symmetric. 

Predicted modal SDCs, natural frequencies, and static deflections at the center of an 
aluminum (Al 4 /i/Al 4 )s and various Gr/epoxy damped plates, (0 4 /i/0 4 )s, ( 0 / 7904)5 and 
(0 2 /90 2 /i/0 2 /90 2 )s are shown in Table 1. All plates have high modal damping values. The 
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laminates with the higher anisotropy in the composite sublaminates, ie. ((yi/0 4 )s and 
(0 4 /i/90 4 )s, exhibit the higher damping in the fundamental mode. The (0Ji/90 4 )s laminate 
exhibits the higher anisotropy variation through-the-thickness among all plates, and this 
seems to contribute to higher modal damping. Interestingly, this laminate has comparable 
static stiffness and natural frequency with the (0 2 /90 2 /i/0 2 /90 2 )s laminate, which illustrates a 
possibility for higher damping values without significantly penalizing other structural 
characteristics. The high damping of the unidirectional composite laminate was attributed 
in part to the lower stiffness of the plate. The composite plates have higher static deflections 
but lower areal density than the aluminum plate, as a result, they exhibit slightly higher 
fundamental frequencies. The damping of the higher modes varies widely and depends on 
the respective mode shape. The variation is more predominant for the composite plates. 

The variations of the first modal damping, natural frequency, and static deflection at the 
center of the plate are plotted in Fig. 2 as functions of the damping layer thickness. The 
thickness of the composite sublaminates was kept constant, therefore, the total laminate 
thickness was varied. The layer thickness t t is normalized by the nominal ply thickness t,= 
0.005 in. The modal SDC is also normalized by the shear SDC of the polymer, and will be 
referred in the next paragraphs , as normalized modal damping. For layer thicknesses in the 
range of practical interest, the laminate with the higher anisotropy through-the-thickness 
(O^i/W^s provides consistently higher damping than the aluminum (Al^i/Al^s and 
(0 2 /90 2 /i/0 2 /90 2 )s laminates. This results are consistent with the expectation that composite 
laminates will induce higher interlaminar shear stresses in the interlaminar layer, in-order 
to balance the elastic anisotropy and nonuniformity of adjacent plies, as opposed to 
traditional isotropic laminates where interlaminar stresses will only balance the bending 
stresses through-the-thickness. For all laminates, the modal damping increased with thickness 
and tends to reach a maximum at higher thickness values. It was observed that for plates of 
lower aspect ratio a/h the peak in the modal damping shifts towards lower thicknesses. 
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Finally, it seems that a critical damping layer thickness exists (in this case approximately 
equal to the ply thickness), such that, damping layers with less or equal thickness will add 
low but significant damping with negligible effects on the structural properties (static 
deflection, natural frequency). 

Fig. 3 shows the predicted variations of modal damping, natural frequency and static 
deflection with the aspect-ratio of the plate (a/h). As was expected, the modal damping of 
all plates increases at lower aspect ratios because shear prevails over flexure. Interestingly, 
the damping reaches a peak at very low aspect ratios and then decreases again. Intuitively, 
at very low aspect ratios through-the-thickness deformations will gradually prevail over shear 
and the effectiveness of the damping layer will be decreased. Nevertheless, the kinematic 
assuptions of the DLDT are not expected to be valid at the regime of very low aspect ratios 
because the present approach assumes negligible through-the-thickness strains. Both natural 
frequency and static stiffness consistently increase with lower aspect ratios. 

Fig. 4 presents the variation of modal damping, natural frequency, and static deflection as 
functions of the FVR of the inner and outer composite sublaminates with respect to the 
damping layers. Fig. 4 demonstrates one unique advantage of composite laminates, that is, 
the additional tailoring capacity they can introduce. In both laminates, the FVR variation of 
the outer sublaminate has a definite effect on the modal damping of the plate. The damping 
increases almost linearly with the FVR, indicating that the most important composite 
parameter is the longitudinal modulus. Since both natural frequency and static stiffness also 
increase with the FVR of the outer sublaminate, the obvious choice is to have outer 
sublaminates of very high FVR. The modal damping is less sensitive to FVR variations of 
the inner sublaminate, and decreases slightly with FVR. 
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Summary 


This paper presented recent research developments on integrated damping mechanics for 
composite laminates with interlaminar damping layers. Damping mechanics for thick 
composite laminates and plates were developed based on a discrete-layer laminate theory. 
The discrete-layer displacement field of variable degrees of freedom enabled the calculation 
of the interlaminar shear strains in the damping layer and composite plies and the 
subsequent prediction of interlaminar shear damping. A semi-analytical method was 
developed for the prediction of the damped dynamic characteristics of composite laminates 
with interlaminar damping layers. 

Applications on Gr/epoxy orthotropic (unidirectional or cross-ply) composite plates with 
layers of damping polymers demonstrated that composite laminates have dynamic 
characteristics comparable to those of geometrically identical aluminum plates. In addition, 
composite plates with strong anisotropy variations between the restraining sublaminates 
exhibited potential for higher damping than geometrically equivalent aluminum plates. The 
effects of thickness aspect-ratio, damping layer thickness and composite fiber volume ratio 
on the modal damping, natural frequency, and static deflection were investigated. Overall, 
the results have demonstrated the potential of embedded interlaminar damping layers in 
composite laminates, moreover, they demonstrated the merits of the developed damping 
mechanics. Future studies will address the effects of damping layers on the laminate strength, 
buckling, and impact resistance. 
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APPENDIX 


Off-axis ply stiffness matrix [E c ] (ref. 13): 

K> = [£J{e c > (Al) 

Transformation matrices: 

m 2 n 2 0 0 0 2 mn 

n 2 m 2 0 0 0 -2 mn 

0 0 1 0 0 0 

0 0 0 m n 0 

0 0 0 -n m 0 

-mn mn 0 0 0 m 2 -n 2 

[rt ]- 1 = [fl(- 0 )] ( A3 ) 

On-axis damping matrix: 

\),J 0 0 0 0 0 

0 0 0 0 0 

0 0 il> /5 0 0 0 

0 0 0 0 0 

0 0 0 0 ^ /5 0 

0 0 0 0 0 i|i w 

The off-axis ply damping matrix is non-diagonal and has the general form: 
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wg = 


4>c// 4 >d 2 0 0 0 y\, cl6 

4* c21 *c22 0 0 0 Xjl c26 

0 0 4; c5J 0 0 0 

0 0 0 *C44 4 * c45 0 

0 0 0 * c5 , 4^ c55 0 

.^c^/ 4* c62 0 0 0 4» c66 


(A5) 
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Table 1. Static and dynamic characteristics of composite plates with interlaminar layers. 
(tj= 0.005 in; numbers in parenthesis indicate mode order) 


(Al 4 /i/Al 4 )s 

(0 4 /i/0 4 )s 

(0/1/904)5 

(0 2 /90 2 /i/0 2 /90 2 )s 

Areal Density, 10' 3 lb/in 2 




7.95 

4.92 

4.92 

4.92 

Static Deflection (Center), in 



0.334 

0.604 

0.496 

0.456 

Modal SDCs 




Mode 1 0.839 (1,1) 

0.978 (1,1) 

0.977 (1,1) 

0.626 (1,1) 

Mode 2 1.656 (1,2) 

0.809 (1,2) 

0.421 (1,2) 

1.703 (1,2) 

Mode 3 2.217 (2,2) 

0.669 (1,3) 

2.771 (2,1) 

1.615 (2,1) 

Natural Frequencies, Hz 




Mode 1 115.7 

119.9 

119.9 

124.0 

Mode 2 265.0 

152.7 

211.8 

272.6 

Mode 3 395.2 

230.1 

354.8 

357.7 
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(a) Typical laminate configuration. 



(b) Kinematic assumptions of discrete layer laminate theory. 
Figure 1 .—Laminates with interlaminar damping layers. 
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Figure 2 —Effect of damping layer thickness on static and dynamic charac- 
teristics of composite plates. Aluminum * (AI 4 /i/AI 4 )s; laminate 1 = 
(0 4 /i/90 4 )s; laminate 2 * ^/QO^i/O^O^s. 
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Figure 3 —Effect of plate aspect rati 
composite plates. Aluminum = (Al 
laminate 2 * (02/902^i/0^902)s, 
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10 

w 


— Aluminum 

— Laminate 1 

— Laminate 2 





,05 .16 .27 .38 .49 .60 


Fiber volume ratio 

Figure 4. — Effect of fiber volume ratio on static and dynamic characteristics of 
composite plates. Laminate 1 = (04/i/90 4 )s; laminate 2 = (0 2 /902/i/02/90 2 )s. 
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